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Abstract 

We derive a second-order differential equation for the two-loop sunrise graph in two dimen- 
sions with arbitrary masses. The differential equation is obtained by viewing the Feynman 
integral as a period of a variation of a mixed Hodge structure, where the variation is with 
respect to the external momentum squared. The fibre is the complement of an elliptic curve. 
From the fact that the first cohomology group of this elliptic curve is two-dimensional we 
obtain a second-order differential equation. This is an improvement compared to the usual 
way of deriving differential equations: Integration-by-parts identities lead only to a coupled 
system of four first-order differential equations. 



1 Introduction 



Precision calculations in high energy particle physics require the computation of quantum correc- 
tions. These can be visualised by Feynman loop diagrams. In this paper we consider a particular 
loop diagram, the two-loop sunrise graph, shown in fig. (OQ). This graph has received in the past 
significant attention in the literature [TTT-IT2TI . Despite this effort, an analytical answer in the gen- 
eral case of unequal masses is not yet known. The state-of-the-art for the two-loop sunrise graph 
can be summarised as follows: In the special case where all three internal masses are equal, a 
second-order differential equation in the external momentum squared and its analytical solution 
are known J9]|. In the general case of unequal masses integration-by-parts identities lfT3l[l4ll 
can be used to relate integrals with different powers of the propagators. In the case of the sunrise 
topology with unequal masses all integrals can be expressed in terms of four master integrals plus 
simpler integrals. This results in a coupled system of four first-order differential equations for the 
four master integrals H). For practical applications this system can be solved numerically fTlfTTll. 

In this paper we reconsider the two-loop sunrise graph in two dimensions with unequal 
masses. We will show that also in the unequal mass case there is a second-order differential 
equation for a single master integral. (The other three master integrals are then given as the 
derivatives with respect to the three internal masses.) This second-order differential equation 
is not obtained from integration-by-parts identities. Setting all masses equal, the second-order 
differential equation agrees with the well-known second-order differential equation for the equal 
mass case. 

How is this second-order differential equation obtained? Our method is interesting in its own 
right and not limited to the special case of the two-loop sunrise diagram. We expect the method to 
be applicable to other loop integrals as well. The starting point is the relation between Feynman 
integrals and periods of motives |[T5l - l2"51l . We view the two-loop sunrise integral as a period of 
a variation of a mixed Hodge structure. The variation is with respect to the external momentum 
squared. In the case of the two-loop sunrise integral we find that the fibre is the complement of 
an elliptic curve. The theory of elliptic curves is well studied. In particular, there is a second- 
order differential equation - the Picard-Fuchs equation - related to a family of elliptic curves. 
The Picard-Fuchs equation expresses the fact that the first cohomology groups of the elliptic 
curves are two-dimensional. The groups are generated by the holomorphic one-form and its first 
derivative with respect to the variation parameter. It follows that the second derivative must be a 
linear combination of these two generators, which leads to a second-order differential equation. 
This is the the sought-after differential equation for the two-loop sunrise graph. 

In this paper we give a detailed account how the second-order differential equation is ob- 
tained. The analytic solution of this equation is beyond the scope of the present paper and will 
be dealt with in a future publication. We expect that a solution can be obtained along the lines of 
refs. [0121. 

This paper is organised as follows: In the next section we define the two-loop sunrise integral 
and recall a few basic facts. Section[3]is devoted to the derivation of the second-order differential 
equation. This section is divided into several subsections. We start with presenting the formalism 
in subsection 13. II The concrete calculations are split into two parts. In the first step we derive 
the coefficients of the homogeneous part of the equation in subsection 13.21 In a subsequent 
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Figure 1 : The two-loop sunrise graph. 

step we obtain the inhomogeneous terms. This is done in subsection 13 .31 For the convenience 
of the reader the complete result is summarised in subsection 13 .41 Finally, our conclusions are 
contained in section @] In an appendix we discuss the relations between the two-loop sunrise 
integrals in D = 2 and D = 4 — 2£ dimensions. 

2 Definition of the two-loop sunrise integral 

The two-loop integral corresponding to the sunrise graph with arbitrary masses is given in D- 
dimensional Minkowski space by 

S(D,p 2 ,m 2 u mj,ml,fj 2 ) = 

J i%% ItS ( K -kl + m 2 ^(-kl + ml)( K -(p-k l -k 2 ) 1 + mf}' 

The corresponding sunrise graph is shown in fig. (QQ)- In eq. (Q} the three internal masses are 
denoted by m\, m 2 and m^. The arbitrary scale /j is introduced to make the integral dimensionless. 
p 2 denotes the momentum squared. This variable plays an important role in our derivation and it 
is convenient to introduce the notation 

t = p 2 . (2) 

Where it is not essential we will suppress the dependence on the masses m,- and the scale /j 
and simply write S(D,t) instead of S(Dj,m 2 ,ni2,m 2 ,fi 2 ). In terms of Feynman parameters the 
two-loop integral is given by 



S(D,t) = r(3-D){f) 3 ~ D f-^to (3) 



a 



with the two Feynman graph polynomials 



11 = X\X2+X2Xi+Xt,X\, 

x\m,\ +X2tn 2 + X3tn 3 ) 11. (4) 
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The differential two-form (0 is given by 

(O = x\d.X2 A dx^ — X2<ix\ A dxj, + x^dx\ A dx 2 . (5) 

The integration is over 

o = {[jq :x 2 :x 3 ] e F 2 \xi > 0,i= 1,2,3}. (6) 

In order to facilitate a comparison with results in the literature we remark that the definition of 
the sunrise integral in eq. ® is in Minkowski space. In a space with Euclidean signature one 
defines the two-loop sunrise integral as 

S'eucl (D,P 2 \m\,rr%,rr%,i?) = 

,^3-d r d D K,d D K 2 1 (7) 

' n% tS (K 2 + m 2 )(K 2 + m 2 )^(P-Ki-K 2 ) 2 + m 2 ^ 

The momenta in Euclidean space are denoted by capital letters, while the ones in Minkowski 
space are denoted by lower case letters. We have the relation 

S &uc i(D,P 2 ,m 2 lJ ml,m 2 i ,//) = S (D, -p 2 ,m\ ,m^mf,// 2 ) . (8) 

Integration-by-parts identities can be used to derive relations between integrals with different 
powers of the propagators 013111411 . Setting 



5 = S(D,t,m\,m\,ml,i? y \ y \,\) = S (D,t,m],n%,ml,i/) , 

51 = S(D,t,m 2 l ,m 2 ,,ml,p 2 ,2, 1,1) = -^ 2 ^5 (D,?,mf,m^,m^ 2 ) 



S 2 = S (D, t, m 2 , n%,n%,i?, 1,2, l) = -a/ 2 ^5 (D,?,m 2 ,m^,m^ 2 ) 



S3 = 5(D,?,m 2 ,m^,m^ 2 , 1, 1,2) = -^ 2 ^|5(D,?,m 2 ,m^,m^,^ 2 ) , (9) 



where 



5(Z), j p 2 ,m 2 ,mi,mf, J u 2 ,Vi,V 2 ,V 3 ) = (10) 
/ 2 N v i+ v 2+v 3 -£» r d D k\ d D k 2 1 

v" > I ~ir~g~~~, 777, ' ~z 77777 ; ; r~; tt^ 



m 3 



is the sunrise integral with arbitrary powers of the propagators, one obtains in this way a set of 
four coupled first-order differential equations for the four quantities So, Si, S 2 and S3. The set 
of first-order differential equations can be found in ref. (TJ. Within the context of integration- 
by-parts identities the four quantities So, Si, S 2 and S3 are referred to as master integrals. In this 
paper we will derive a single second-order differential equation for the quantity So- 
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It will be convenient to derive this equation in D = 2 dimensions. Dimensional recurrence 
relations [[27H301 can be used to obtain the result in D = 4 — 2e dimensions. Since the application 
of the dimensional recurrence relations involves a subtlety, we have addressed this issue in an 
appendix. Working in D = 2 dimensions has two advantages: First of all, the integral S(2,t) 
is finite and does not require regularisation. The second advantage is given by the fact, that in 
D — 2 dimensions eq. © reduces to 



5(2,0 = ffj- (ID 



In D — 2 dimensions the Feynman integral depends only on the second Symanzik polynomial jF, 
which occurs in the denominator, but not explicitly on the first Symanzik polynomial U. 

Let us briefly discuss the simpler case of equal masses. Setting mi = mi — m?> = m it has been 
shown by Laporta and Remiddi that the system of four coupled first-order differential equa- 
tions reduces to a single second-order differential equation. In D = 2 dimensions this equation 
reads: 



d 2 (3t 2 -20tm 2 + 9m 4 ) d t-3m 2 



dt 2 t (t - m 2 ) (t - 9m 2 ) dt t(t-m 2 )(t-9m 2 ) 



5(2,0 = -6 " 



t(t-m 2 )(t-9m 2 )' 

(12) 



In this paper we will show that this second-order differential equation is not an artifact of the 
special case m\ = m% = = m, but has its origins in the interpretation of the zero set of 'J 
within algebraic geometry. 



3 Derivation of the second-order differential equation 

In this section we show that the two-loop sunrise graph with arbitrary masses has a second order 
differential equation similar to the one in eq. ([12)) in the case of equal masses. In this section we 
assume that all masses are positive and rational. 

3.1 Formalism 

Our starting point is eq. (fTTI) : 

5(2,0 = f?jj = J (O t , (13) 
o a 

where we have set (D t = /j 2 (o/^T. We interpret this Feynman integral as a period of a variation 
of a mixed Hodge structure (VMHS), varying with t. Although the fibres depend on the masses, 
the rank of the VMHS is independent of them and so the differential equations should be of the 
same complexity in the unequal mass case compared to the equal mass case. We denote by X 
the set of points ([jci : X2 : xj] ,t) & P x A*, for which = 0. Here, A* is an open subset of C. 
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We denote the fibre over t by X t . The second Symanzik polynomial [F is of degree three in the 
variables x\, X2 and X3 and the zero set defines a family of elliptic curve in P 2 depending on t. 
Then Q) t G H 2 (F 2 \X t ) for all t, but o i/2(P 2 \Xr)> i- e - is not a cycle, for two reasons. Firstly, o 
intersects X t and secondly o has a boundary. The second fact being obvious, we now deal with 
the first. 

Lemma 1 Let to := {m\ +ni2 + and C <to be the complex numbers with the line {x G R | 
x>to} removed. For any t G C <fo ?/ze c/zam of integration o intersects the graph hypersurface 
X t precisely in the three points [1:0:0], [0:1:0] and [0:0: 1]. 

Proof 1 We have f = — t X1X2X3 + {m\x\ +m\x2-Vm\x-s){x\X2 J cX\x? > + X2X3). First observe that 
the boundary of a intersects X t precisely in the three points stated. We have to show that the 
intersection of X t with the inner points of a is the empty set for t G C< fo . This is obvious for 
t G C\ M. Now let t = (mi +m2 +m3) 2 — 8, with 8 G M>o- We restrict to the affine open x\ = 1 
and obtain the function 

*F = —((mi +m2 + mT,) 2 — 8)^2^3 + (m\ -\-m\x2 J rm\xi)(x2 +X3 +X2X3). 

We have to show that the equation 

(mi+m2 + m3) — 8 = (m\ +m 2 X2 + m 3 X3)( 1 hi) 

X2 X3 

has no positive real solution. Now 9(^2,^3) := (m\ + m\x2 +m 2 X2,)(-^ + ^- + 1) is a continuous 
function from U := IR>o x R>o to M>o which tends to infinity, when X2 or xj tend to zero or 
infinity. Hence the set {x G U \ q>(x) < C} =: K C U is compact and cp has it's global minimum 
on K. We easily find the global minimum to be unique, namely the point (^2,^3) = ^-). 

Now we have <p(^, ^) = («^i +m2 + m3) 2 which proves the Lemma (the other two affine opens 
give the same solution). □ 

In the following we will assume t G C< ?0 . The differential equation which we derive will be 
valid in the region C <tQ . Note that p 2 = to = (m\ +m2 + mj,) 2 is the physical threshold. The 
two-loop sunrise integral for values of p 2 above the threshold can be obtained from the solution 
of the differential equation by analytic continuation with the help of Feynman's z'£ prescription. 

Now let P — > P 2 be the blow up of P 2 in the three points of Lemma [IJ We denote the strict 
transform of X t by Y t and the strict transform of o again by a. In the particular example of the 
two-loop sunrise graph we are in the lucky situation that X t is isomorphic to Y t for generic t - 
both are smooth elliptic curves. In P we have oHF, =0. Now let Bq := {X1X2X3 = 0} C P 2 and 
B it's total transform. Clearly the boundary of o is contained in B. We now find 

H t :=H 2 (P\Y t , B\BC]Y t ) 

to be the right mixed Hodge structure, i.e. (0, G H 2 (P \ Y,, B \ B n Y t ) and o G H 2 (P \ Y t , B \ B n Y t ) . 
Note that this is very similar to the work of Bloch, Esnault and Kreimer [15]. The convergent 
Feynman-Integral S(2,t) is a period of H 2 (P \Y t , B \ BnY t ). We compute it's Picard-Fuchs 
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equation. 



In the following we will denote a generic fibre by X, resp. Y , dropping the subscript t. Since 
Y is smooth we can apply the Gysin sequence which is a sequence of mixed Hodge structures 
and reads 

— >H\P\Y) — >H°(Y) — >H 2 (P) — >H 2 (P\Y) — >H l (Y) — > 0, 

where we have used H l (P) = H 3 (P) = 0. Indeed we have H k (P) = H k (P 2 ) foxk^l and 
H 2 (P) = 7t*( [w] ) © ZEi © ZE 2 © ZE 3 , where w is a generator of H 2 (P 2 ) and the Ej correspond to 
the exceptional divisors of the blowup (which are P 1 *). But i* maps H°(Y) isomorphically onto 
Ti*(w), such that we get H l (P\ Y) = and a short exact sequence 

> ZEi®ZE 2 ®ZE 3 > H 2 (P\Y) H l (Y) > 0. (14) 

This sequence is split as a sequence of mixed Hodge structures via 

H 2 (P\Y) H\Y) 
H 2 (F 2 \X) H l (X). 



The elliptic curve X has a unique holomorphic one-form (unique up to exact forms), which 
reads in Weierstrass normal coordinates y if we restrict to the affine open z = 1. The Picard- 

Fuchs equation of y 6 H l (X) can easily be computed. We find the Picard-Fuchs operator 

i,2, = £ +<l(!) s +6W - (16) 

with rational functions a(t) and £>(?). 

This is also the Picard-Fuchs operator of the Feynman form oo f E H 2 (P\ Y) due to the splitting 
of sequence ([141) and the flatness of the system 1E\ © 7LE2 © 7LE 3 . So for any cycle ^ in H2(P\ Y) 
we have 

z/ 2 )(|o),) = o. 

Since the domain of integration o is not a cycle in #2 (P\Y), as explained above, we now 
pass to the relative setting. There is the long exact sequence of relative cohomology 

— >H\B\BnY) — >H 2 (PXY, B\BHY) — >H 2 (PXY) — >H 2 (B\BnY), 
where we have used H l (P\Y) = 0, as shown above. But now we find 

BnY = {[0: : 1],[-^ : 0: 1],[1 : : 0]} U { Pl ,P2,P3}, 
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Figure 2: B\ B D Y. The lines are P s. The circles indicate, that one point is removed from each 
copy of P 1 . 



where pi is a point on the exceptional divisor 

For BXB fl y we get the picture of figure [2] Using the Mayer-Vietoris sequence we get 

H°(BXBnY) = Z, 

fl^BVBny) = Z(-1), and 

H k (BXBnY) = 0,for£^ 1,2. 

Summarising, we get the short exact sequence 

o > z(-i) > // 2 (p\y, 5\5ny) > h 2 (pxy) ► o (17) 

of mixed Hodge structures. 

Proposition 2 The Picard-Fuchs operator of the Feynman integral S(2,t) is 
L := — ^Mz/ 2 ), where g(t) is a rational function oft with complex coefficients. 

Proof 2 We have seen some one-form pY This amounts to the identity 

L^(J(D t ) = fdb = j % =:g(t). 

a a da 

From sequence ( fT7|) we know that the Picard Fuchs operator of S(2,t) can have order three at 
most. Obviously L(S(2, t)) = 0. It remains to show that L^> (S(2,t))^ 0. Our explicit calculation 
shows that this is indeed the case. □ 
From 



U 

we obtain immediately 

dt 2 ' " X ' J dt 



fdh) = (18) 



-^+a(t) — + b(t) 



S(2,t) = g(t). (19) 



This is the sought-after second-order differential equation. It remains to determine the coeffi- 
cients a(t) and b(t), as well as the inhomogeneous term g(t). 
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3.2 Calculation of the homogeneous part 



In this subsection we outline the calculation of the coefficients a(t) and bit) in the homogeneous 
part of the differential equation. From eq. (fT5l) it follows that it is sufficient to consider H (X). 
We recall that the variety X is defined for fixed t as the zero set in P 2 of the second Symanzik 
polynomial 'J: 



— XiX2X3t+[xirnl+X2m2+X^mfj(xiX2+X2X3+X3Xi) = 0. 



(20) 



This polynomial is of degree 3 in the variables x\ , X2 and X3 and defines an elliptic curve. By an 
birational change of coordinates the defining equation can be brought into the Weierstrass normal 
form 



y 2 z - x 3 - a 2 (t)xz 2 - a 3 (t)z 3 
In the chart z = 1 this reduces to 

y 1 — x 3 — a2(t)x — a?,{t) = 
In this coordinates H l (X) is generated by 

dx d 
ri = — and r) = — ri 

dt 



0. 



y 



(21) 



(22) 



(23) 



Since H (X) is two-dimensional it follows that f\ 
f|. In other words we must have a relation of the form 

f\ + a(t)r\ + b(t)r[ = 0. 



d T| must be a linear combination of T| and 



(24) 



It is convenient to bring this equation onto a common denominator. Doing so and carrying out 
the derivatives with respect to t we have 

/ q \ 2 dx 

T| = [X +a2X + d3) -j, 



1 dx 

-- (d2X + d3) (x 3 + a2X + a-i) — , 
L y 

1 3 

— - (azx+di) (x 3 +a2X + as) + - (d2X + ds) 



dx 

7' 



(25) 



The numerator of eq. (1241) is then a polynomial of degree 6 in the single variable x. Since we 
work in H l (X), we can simplify the expression by adding an exact form 



ji-l 



n — - x + In — - CI2X + n«3 



dx 

7' 



(26) 



This allows us to reduce the numerator polynomial from degree six to a linear polynomial. The 
two coefficients of this linear polynomial have to vanish, on account of eq. (l24l) . We obtain 
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therefore two equations for the two unknown parameters a(t) and b(t). Solving for a(t) and b(t) 
we find 



a(t) 



Pi(t) 
Po(t)' 



b(t) 



P2(t) 

Po(t)' 



(27) 



with 

Pi(0 



Pi{t) 



Po[t) 



9t 6 - 32M m t 5 + (37M 200 + 10M X 10 ) t A - (8M 300 + 56M 210 + 144M i n) ? 3 

- ( 13M 400 - 36M 3 io + 46M 220 - 124M 2 n ) f 2 

- (-8M500 + 24M 4 io - I6M320 - 96M 3 n + 144M 22 i) t 

- (M 600 - 6M 5W + 15M 420 - 2OM330 + I8M411 - 12M 32 i - 6M222) , 
3t 5 - 7M l00 t 4 + (2M 200 + 16Mno) t 3 + (6M300 - 14M 2 io) t 2 

- (5M 40 o - 8M310 + 6M220 - 8M211 ) t + (M500 - 3M 4 io + 2M 320 + 8M311 - IOM221) 

,2" 



2 2 

t — (mi+m2 + m^) t — (— m\ + ni2 + m^) 



(mi — ni2 + tn3) 



t — (m\ +ni2 — m^) 



[3t 2 -2M m t-M 2 oo + 2M no ] . 



(28) 



In order to present the result in a compact form we have introduced the monomial symmetric 
polynomials Mj^juju m me variables m 2 , and m 2 . These are defined by 



(29) 



where the sum is over all distinct permutations of (A4 , A, 2 , A3). A few examples are 

M100 = mf + mj + m 2 , 
Mm = mim 2 m3, 



M 2 
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42 42 42 42 42 42 

m l m 2 + m 2 m 3 + m 3 m i +m 2 mi +m 3 m 2 + mim 3 . 



(30) 



3.3 Calculation of the inhomogeneous part 

In this subsection we compute the inhomogeneous part g(t). From proposition 2 we first seek a 
one-form p\, such that 

d& = L^((D t ) (31) 

We make the ansatz iBTfl 

Pf = — 7^-ppj[(x2q3-X3q2)dxi + (x3qi-xiq3)dx2 + (xiq2-X2qi)dx3], (32) 

where gi, q2 and ^3 are polynomials of degree 4 in the variables x\, %2 and X3. The most general 
form is 

Gi = 

JO v.4,J0 J , J) 4 1 (0 v 3„ 1 JO V 3 V , 1 JO VlV 3 , JO Y 3 , JO V 3 V , JO „ v 3 
c 400 1 ' 040 2 ~r L 004 A 3 ~r L 310 A 1 A 2 ~r c 301-*l-*3 1" C 130 A 1 A 2 ~r L 103 A 1 A 3 C 031 A 2 A 3 1" c 013 2A 3 

4- JO v 2 v„v, 4- JO v,v 2 vo-uJ0 v ,v.v 2 -uJ0 v 2 v 2 4-J0 y 2 ^ 2 , JO v 2 v 2 (W\ 

"^ 21 1 1 2-^3 ~r t- 121"^1 2-^3 "T t- j ^2-*l-*2-*3 ~r ^220 1 2 ~r ^202 1 3 022 2-^3 • Ip-V 
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We would like p\ to be finite on the boundary do. This implies 

c (l) _ c (l) _ c (2) _ (2) _ (3) _ (3) _ 

c 040 ~~ c 004 — c 400 ~~ c 004 ~~ c 400 ~~ c 040 ~~ w " \ DH > 

The remaining 39 coefficients c% are found by solving the linear system of equations obtained 
from inserting the ansatz into eq. (TJTI) . The solution of this linear system is not unique, cor- 
responding to the fact that p\ can be changed by a closed one-form. The solutions for the co- 
efficients are rather lengthy and not listed here. In the next step we integrate p\ along the 
boundary do to get g(t): 

(35) 



3a 

2 



Note that the integration is in the blow-up P of P . We obtain 



fA P3 

git) = —7^, (36) 
P0{t) 

with 

P3 (t) = -m 4 + 24M m t 3 + (4M 200 - 40M no ) t 2 + (-SM 3 oo + SM 2W + 4SM m )t 

2 

+ ( -2M 400 + 8M 3 io - 1 2M 220 - 8M211 )+2c(t,m u m 2 ,m 3 )\n-t 



2 2 
+2c (t,m,2, m3, mi ) In + 2c (t, m 3 , mi , m 2 ) In (37) 



and 



c(f,mi,m 2 ,m 3 ) = 

(-2m? + m 2 + m 3 ) ? 3 + (6m{ - 3m 2 - 3m 3 - lm\m\ - lm\mk, + 14m 2 m 3 ) t 2 

I f\ fit f\ AO AO OA OA A A\ 

+ ( — 6mj + 3m 2 + 3m 3 + 1 lmim 2 + 1 lmim 3 — 8mim 2 — 8mim 3 — 3m 2 m 3 — 3m 2 m 3 J t 
+ f 2m? -m\-m\- 5m\m\ - 5m\ml + m]m\ + mjm^ + 4m 2 m 3 + 4m 2 m 3 
+3mX + 3«X - 6m\m\ + 2m?m 2 : m 2 , - m\m\m\ - m\m\rr§ . (38) 
The coefficients c(t, m ; , rrij, nik) of the logarithms of the masses vanish for equal masses. 

3.4 Final result 

In this subsection we summarise our results. The two-loop sunrise integral S(2j) with unequal 
masses satisfies a second-order differential equation. This second order differential equation is 
given by 



d 2 pi(t) d p 2 (t)' 
dt 2 po(t) dt po(t) 
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5(2,0 = (39) 



where po(t), p\{t), pi{t) and p^it) are polynomials in t. The polynomials po(t), p\{t) and pi{t) 
are defined in eq. (|28l) . while P3(t) is given in eq. (1371) . In the special case of equal masses 
m\ = m,2 = m-s = m, eq. (l39l) reduces to eq. dT2]) . 

4 Conclusions 

In this paper we have shown that the two-loop sunrise integral in two dimensions with three 
arbitrary masses has a second-order differential equation. This differential equation is derived 
in a new way: We view the Feynman integral as a period of a variation of a mixed Hodge 
structure, where the variation is with respect to the external momentum squared. We expect 
our technique to be applicable to other Feynman integrals as well, in the following sense: In 
seeking a differential equation with respect to an external invariant for an unknown Feynman 
integral we view the Feynman integral as a period of a variation of a mixed Hodge structure. 
The fibre will be the complement of an algebraic variety. The algebraic variety is given by the 
zero set of the Feynman polynomial J and depends therefore on the specific Feynman integral 
under consideration. In the specific case of the two-loop sunrise integral the algebraic variety is 
an elliptic curve, while for other integrals of interest we expect to encounter different algebraic 
varieties. 
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A Dimensional recurrence relations 

Although the focus of this article lies on the two-loop sunrise integral in two dimensions, we 
would like to give in this appendix some indications how the integral in two dimensions can be 
related to the one in D = 4 — 2e dimensions. In the unequal mass case there are some subtleties, 
which are discussed in this appendix. 

Dimensional recurrence relations can be used to relate the two-loop sunrise graph in D = 
4 — 2e dimensions to the result in D = 2 — 2e dimensions. Let us introduce an operator i + , which 
raises the power of the propagator i by one, e.g. 

l + 5(D,?,mf,w^w^Vi,v 2 ,v 3 ,M 2 ) = 5(D,?,m?,m^m^Vi + l,v 2 ,v 3 ,^ 2 )- (40) 
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The integral S(Z),£,mi,m2,m3,Vi,V 2 ,V 3 ,/r) is defined by 

S (D,p 2 ,m 2 l ,ml,ml,v l ,V 2 ,V3,V 2 ) = (41) 
/ ?nv-d f d D k\ d D k2 

= (a< 2 ) 1 



in? in? (_ k 2 +m iyr (_ k 2 +m 2yi ^_( p _ h _ k2) 2 + 

r( Vl )r(v 2 )r(v 3 ) ^ j 7 Xl * 2 * 3 w ' 



V3 



where we have set v = Vi + V2 + V3. 

The starting point for the dimensional recurrence relation is given by H27H28H 

5(D-2,f,m?,m^m^Vi,V 2 ,V 3 ,A/ 2 ) = U (vil + ,V 2 2+,v 3 3 + ) 5(D,/',m 2 ,m^m^Vi,V2,V3,^ 2 ). 

(42) 

Eq. (1421) is easily verified by writing out the right-hand side: 

U (vil + ,V 2 2 + ,v 3 3 + ) 5(A?,m 2 ,m^m^Vi,V 2 ,V 3 ,A/ 2 ) = 
= ViV 2 S(D,?,m 2 ,m 2 ,m 2 ,Vi + 1,V 2 + 1,V 3 ,^ 2 ) +V 2 V 3 5(D,?,m 2 ,m 2 ,m 2 ,Vi,V 2 + 1,V 3 + 1,/i 2 ) 
+V 3 Vi5(D,?,m 2 ,m^m^Vi + 1,V 2 ,V 3 + 1,/i 2 ). (43) 

Using eq. (|4TT) we can express 

ViV 2 5(D,?,m 2 ,m^m^Vi + 1, V 2 + 1, V 3 „u 2 ) = 



r(vi)r(v 2 )r(v 3 ) v ^ y J 1 2 3 !F v + 2 - D 



and similar for the two other integrals on the right-hand side of eq. (1431) . Adding up the three 
integrals we find 



1i(vil + ,V 2 2 + ,v 3 3 + )5(A^mf,m|,m^,Vi,V 2 ,V 3 ,A/ 2 ) = 

T(v + 2-D) , 2 ,v+2-D f Vl _! V2 _! v,_i ^ +3 -§ Z) 



(u 2 Y +2 - D fr^-K^-K^- 1 U 2 

\H J J x l x 2 x 3 yrv+2-D 



(0, (45) 



r(vi)r(v 2 )r(v 3 ) v " y y if 

which is just S(D — 2,f,m 2 ,m2,m 2 ,Vi,V 2 ,V 3 , J u 2 ). 

Integration-by-parts identities |[T3l[T4l[32l are of the form 

, 2 vv-d r<Phfhd_ ^ 

^ j 7 m§ ftcf dtf (-^ 2 + m 2 ) Vl (-* 2 + mf) V2 (-fc 3 2 + m2) V3 U ' 1 0j 

with G {&i,fc 2 }, <? G {p,&i,£ 2 } and fc 3 = p — &i — fc 2 . Integration-by-parts identities can be used 
to derive relations between integrals with the same topology but different powers of the propaga- 
tors and integrals with simpler topologies. Using integration-by-parts identities we express the 
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right hand side of eq. (|42l) in terms of the master integrals So, Si, S2, S3 and simpler integrals, 
which are products of one-loop tadpole integrals. The one-loop tapole integral is given by 

/ ? y\ / ?n1-# rd D k 1 / D\ (m 2 \^~ X 

r (D ,mV) = (, 2 ) ! /^(=FT^)= r ( 1 -2)UJ ' (47) 

Inverting a linear system of equations we can express the D dimensional integrals in terms of the 
(D — 2) dimensional integrals. Since the involved expressions are rather lengthy, this is usually 
done by computer algebra. Specialising to D = 4 — 2e we obtain the relation 

S (4 - 2e, = c (e, ?) S (2 - 2e, ?) + c\ (e, t) Si (2 - 2e, ?) + c 2 (e, t) S 2 (2 - 2e, 

+c 3 (£,0S 3 (2-2e,0 + r(£,0- (48) 

The dependence of the coefficients cq, c\, c 2 , C3 and r on the masses mi, m 2 and m3 is not shown 
explicitly. All quantities in eq. (|48l) can be viewed as a Laurent series in £. The Laurent series of 
S(4 — 2e, t) starts at l/£ 2 and one is usually interested in the pole terms and in the £°-term. The 
quantities S(2 — 2e,t), Si (2 — 2e,t), S 2 (2 — 2e,f) and S3 (2 — 2e,f) are finite, and their Laurent 
series start at 8 . If it is the case that also the Laurent series of the coefficients cq, c\, C2 and C3 
start at £°, eq. (1481) would reduce to 

S(4-2£,f) = co(0,OS(2,0+ci(0,OSi(2,0 + c 2 (0,OS 2 (2,0 + C3(0,OS3(2,0 

+r(E,0 + O(e), (49) 

and it would be sufficient to know S(2,t), S\(2,t), S 2 (2,?) and $3(2,?) in order to determine the 
£°-term of S(4 — 2e, t) . The ultraviolet poles of S(4 — 2e, t) would then be entirely given by the 
tadpole contributions contained in the coefficient r. Unfortunately it turns out that for this choice 
of master integrals the coefficients cq, c\, c 2 and C3 contain spurious poles in £. In order to avoid 
the computation of the E 1 -terms of the master integrals it is advantageous to choose a different 
basis for the master integrals. A basis which does not lead to spurious poles is given by 

h(D,t) = S(D,t,mim 2 2 ,mlp 2 ), (50) 

I 2 (D,t) = v 2 ^rh(D,t), 
at 

r fn A ( 2\2-d f d D k l d D k 2 2p-k\ 
h{D,t) = [fi ) 



fn§ ^kj + m 2 ) (-kl + m 2 ) ^-(p-k l -k 2 ) 2 + ml 



h(D,t) = (p 2 ) 



2,2-D [d D ki d D k 2 2p-k 2 



i%2 ini (_^2 +m 2^_^2 +m 2^_^_^_^2 + m 2 



The second-order differential equation can be used to determine h{2,t) and 7 2 (2,?). In addition, 
the integrals h(2,t) and h(2,t) are required. These are given by 

4(2,0 = V/^, 

/ 4 (2, ( ) = V/^- (51) 
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Only one integral of these two needs to be known, the other follows then from the substitution 
m\ -H- m,2. We can also derive a differential equation for 1$ (or I4) with the methods of this article. 
The actual calculations will be slightly more complicated due to the extra factor X2X3 / 11 (or 
X1X3/U). The integrals I3 and I4 will be treated in detail in a future publication. 
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